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1 Introduction 



In our previous paper [T] we have shown that HERA F2 data, at low x, can be very well 
described by the gluon density constructed from the discrete spectrum of eigenfunctions of 
the BFKL kernel, i.e. from the pomeron wave functions. This first successful confrontation 
of the BFKL formalism [5] with data led to the unexpected question as to whether the HERA 
data are sensitive to the Beyond Standard Model (BSM) effects. These effects, although only 
present at scales that are much higher than the region of HERA data, can nevertheless affect 
the quality of the fits to data because the shape of many of the contributing eigenfunctions 
has an apparent sensitivity to the BSM effects. This apparent sensitivity is due to the fact 
that the support of eigenfunctions extends to very high transverse momenta where BSM 
effects have to be present. Since the eigenfunctions are constructed in a global way, i.e. 
the behaviour of the eigenfunctions at energies way above the threshold feeds into their 
behaviour at low energies, these eigenfunctions will be sensitive to any BSM physics. 

In this paper we investigate whether this possible sensitivity to BSM effects indeed exist, 
as it also depends on an adequate treatment of the infrared boundary condition. As a 
popular example of BSM effects we have chosen the N=l supersymmetry. For this purpose we 
modified the beta function and the kernel of the BFKL equation to include the contributions 
from the superpartners and confronted the modified gluon density with data. 

The paper is organized as follows; in Section 2 we give a brief summary of the proper- 
ties of the discrete pomeron solution to the BFKL equation. In Section 3 we describe the 
construction of the infrared boundary condition and the changes introduced by the two-loop 
running of the coupling. In Section 4 we discuss the effects of the supersymmetric changes 
of the /3-f unction and of the eigenfunctions of the BFKL kernel. In Section 5 we present and 
discuss the results and in Section 6 we describe the properties of the determined infrared 
boundary. In Section 7 we give a summary. 

2 The Discrete Pomeron Solution to the BFKL Equa- 
tion 

In this section we give a brief summary of the properties of the discrete pomeron solution to 
the BFKL equation, described in detail in [1]. 

The BFKL amplitude for the scattering of high-energy gluons with transverse momenta 
k and k', is a Green function constructed from the discrete eigenfunctions of the BFKL 
kernel, i.e. the solutions fn{k) to the equation 



where as(A;^) is the strong coupling which runs with the magnitude, k, of the transverse 
momentum of one of the gluons. This running leads to an oscillatory eigenfunction, /„, whose 
frequency, in the semiclassical approximation, depends on the transverse momentum, /c. 




(2.1) 
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so that the eigenvalues, a;„, are given in terms of the LO and NLO characteristic functions 
of the oscillation frequency u, 

^ = I jXol^^j+l 1 Xi{'^) + ---, (2.2) 

where for the moment we have ignored the resummation of coUinear divergences in the NLO 
characteristic function [5]. The frequency depends on u and decreases as k increases, reaching 
a critical point kcru at I'lkcru) = where it changes from real to imaginary values. Below 
kcrit the eigenfunction has an oscillatory behaviour but above kcrit it decreases exponentially 
with In k. The matching of the phases immediately below and above the critical point fixes 
the phase of the oscillations at kcrit to be — 7r/4 (this being the phase of the Airy function, 
which was shown in ref. [HEllT] to provide a very good approximation to the eigenfunctions) . 
This phase, 4>{k), at any lower value of k is then determined by integrating the fc-dependent 
frequency from kcrit to k, namely 

(P{k) = -- + 2 [ ^"\{k)d\n{k). (2.3) 
4 .1 u 



For a given value of oo the phase, rj vr, of the oscillation at some infrared transverse 
momentum, ko, is determined from the perturbative BFKL dynamics (with running cou- 
pling), through eqs.f l2.1|2.2p . We make a very general assumption that the infrared (non- 
perturbative) properties of QCD impose some phase at k^, defined up to an ambiguity of 
nn, which can also depend on u. We find then that we can only match this phase to that 
determined from eq. fl2.3p for one value of u, for each integer n, where n corresponds to 
the number of turning points of the eigenfunctions. This leads to the quantization of the 
spectrum (i.e. discrete pomeron poles), in keeping with the predictions of Regge theory. 

Before a comparison can be made with the measured structure function, F2, it is necessary 
to convolute this Green function with the impact factors for the virtual photon at one end 
and for the proton at the other, (see Section 6 of ref. pQ). The impact factor for the virtual 
photon is calculable in perturbative QCD and has support, which is peaked at transverse 
momenta of the order of the argument of the structure function, whereas the proton 
impact factor cannot be so calculated and is assumed to have a simple form with support 
up to C(l) GeV. 

One of the main results of ref. pQ was that a very good quality fit to HERA-F2 data [2] 
(with above 8 GeV^) is obtained taking a very simple form for the dependence of the 
above-mentioned (non-perturbative) infrared phase, rjn, on eigenfunction index n. It was 
found that in order to obtain this good description of data, it was necessary to take around 
120 eigenfunctions of the BFKL kernel. 

Although the oscillation frequency varies with transverse momentum, k, the period. A, 
of oscillation (in ln{k)) defined by 

/•ln(fc)+A 

2 / iy{k')d\n{k') = 27r (2.4) 

iln(fc) 
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Figure 1: The running of across a threshold for N=l SUSY at 10 TeV 

turns out to be roughly constant (A ~ 8), beyond the first two turning points. This means 
that the main difference between the n*^ and (n + 1)*'* (for n > 2) eigenfunction is that the 
latter has one more half period, which leads to a rapid increase in the critical momentum, 
kcrit with eigenfunction index n 

k,„t ~ c-e^", (2.5) 

where c is a constant of the order of Aqcd- For the first eigenfunction the value of kcrit is 
0(10 GeV). It therefore follows that kcrit rapidly exceeds the threshold for most postulated 
theories beyond the Standard Model. On the other hand, if a threshold for new physics 
does indeed exist, the oscillation frequency is affected above this threshold and consequently 
the oscillation phase at all lower transverse momenta will be altered thereby affecting the 
matching of the phase to the phase imposed by the infrared dynamics of QCD. This in turn 
modifies the pomeron spectrum, a;„. It is in this sense that a modification of the high-energy 
behaviour of the eigenfunctions "feeds" into the low-energy behaviour. 

This immediately posed the question as to what the effects would be on the quality of 
the fit, if there were some new physics far above the energy scale of HERA. 

3 The Infrared Boundary 

In ref. [1], we defined the infrared boundary as a phase condition at the lowest possible 
value of the transverse momentum, k = k^, which can be safely reached by the perturbative 
calculation. To make this value as close as possible to Aqcd we considered only the one- loop 
running of the coupling. This gave a value of k^ = 0.3 GeV, which corresponds to ~ 0.7. 
The reason for running the coupling at one-loop only was that in principle this is the same 
order of perturbation theory as the NLO characteristic function, xi 

However, given that we modify eq. fl2.ip by resumming all the large corrections in xi 
using the technique of ref. [S], it is more appropriate to take the /3-function to two- loop order 
which is what we use in this paper. 
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When we do this, we are faced with a problem - namely that we cannot run the coupling 
below an "infrared" scale fco = 0.6 GeV, which corresponds to ~ 0.7 (at the two loop 
level), without approaching the Landau pole too closely. On the other hand, the infrared 
boundary conditions are to be imposed at a transverse momentum of order Aqqd- Moreover 
we need to know the eigenfunctions below ko in order to perform a convolution with the 
proton impact factor, which has support mainly below the ko value. Therefore, guided by 
the behaviour of the eigenfunctions in perturbative QCD, we continue them down to a lower 
momentum k^, which should be of order Aqcd, using the extrapolation of the phase 0„(A;) 

Mk) = Mko)-2i^'jn(i-y (3.1) 

where for each eigenfunction, with index n, z/° is the frequency of the oscillations near 
k = ko [1]. We have assumed that this frequency is constant below ko, an assumption which 
is correct for sufficiently small ko, at least for the leading order BFKL kernel (see [7J). Any 
deviation from constant frequency should have a negligible effect as we are only extrapolating 
over a small range in gluon transverse momentum. The numerical values of z/° are obtained 
by inverting the eigenvalue equation ( 12. 2p . modified according to [5]. 



4 N=l Supersymmetry at various Thresholds 

The "new physics" that we investigate in this paper is the popular = 1 supersymmetric 
extension of the Standard model above a given threshold fcy, which for simplicity we assume 
to be a common mass threshold for all superpartners. Below this threshold the running of 
the coupling is governed by the /3-function to two-loop order 

where for the case of QCD, Ca = 3, Cp = 4/3 and is the number of active flavours. 
Above the threshold, the beta function is given by 

/3> = -g {3Ca - -/) - {qCI + (-^ + 2C,^ n,) . (4.2) 

This leads to a "kink" (discontinuity in the derivative) in the running of as at the threshold 
for N=l SUSY, which can be seen in FigUl 

The fact that the coupling runs more slowly above the SUSY threshold means that 
the oscillation frequency varies more slowly with k and therefore the critical transverse 
momentum, kcrit (where u = 0), is pushed out further away. Thus, for example, in the 
case of a SUSY threshold at 10 TeV, if we assume that the phase of the oscillations is 
zero at ko = 0.6 GeVH the first two eigenfunctions are identical as the critical momentum 

^This phase is merely an example designed to illustrate the difference in the behaviour of the eigenfunctions 
for the Standard model and the N—1 SUSY model. In practice these phases are determined from the fit to 
data and are different for the two models. 
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Figure 2: Oscillation frequencies as a function of gluon transverse momentum for various 
eigenfunctions. The left-hand pane is the case of the Standard Model and the right-hand pane 
is the case of N=l SUSY above a threshold of 10 TeV. For the purpose of this comparison 
it has been assumed that the infrared phases are the same in both cases. 



is below the threshold, whereas for the third eigenfunction the critical momentum is at 
kc„t = 1-2 X 10^ GeV in the case of the SM but at k^ru = 1-3 x 10^ GeV in the case of SUSY. 
This can be seen from Figj2l 

Furthermore the NLO characteristic function, Xi{'^) acquires an additional contribution 
19J of 

/ N TT^ sinhfTTZ/) /II 9\ , , 

- 32 .(1 + .-)cosh-(..) (t + ^-V 
from the octet of Majorana fermions (gluinos), and 

/ N vr^ rif sinhfvrz/) /5 r,\ ^, 

from the squarks. Note that it is this discontinuity in xi at the SUSY threshold which 
is responsible for the discontinuities in the frequencies at threshold and not the change in 
the rate of running of the coupling, which remains a continuous functioiJl. The change in 
frequency thus compensates for the change in the characteristic function in order to ensure 
that the eigenvalues, cOn remain unchanged as one passes through the threshold B 

The contribution, 6xi, of these additional terms to xi is shown as a function of frequency 
in Fig. [3] where it can be seen that this is a rapidly decreasing function, which explains why 
the discontinuities in frequency at threshold are much larger for the lower eigenfunctions for 
which the frequency at threshold is lower. 



^ A similar smaller discontinuity can be seen at around 3 GeV. This corresponds to the c-quark threshold. 
There are analogous, even smaller, discontinuities at the b-quark and t-quarks thresholds 

■^The discontinuous changes in frequency are due to the fact that the change in characteristic function is 
imposed at a threshold in its entirety. A determination of the NLO characteristic function which accounted 
for the mass of internal particles would smooth out these discontinuities. 
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Figure 3: Increase in the NLO characteristic function, xi as a function of frequency v 



0.3 



0.25 



0.2 



0.15 



0.1 



0.05 




2 4 6 8 10 12 14 16 18 20 

n 



Figure 4: The first 20 eigenvalues in the case of the Standard model (triangles) and SUSY 
at a threshold of 10 TeV (squares) 
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The differences in frequencies also affect the magnitude of the eigenvalues, beyond the 
first two, as can be seen in Fig. HI The simplest way to understand this is to consider 
the value of at kcru for the two models. Although kcru is an order of magnitude larger 
for the SUSY model, the fact that the coupling runs more slowly actually means that the 
value of the coupling at kcru is slightly larger for the SUSY model. Moreover the NLO 
characteristic function is larger in the SUSY model. These two effects combine to produce 
somewhat larger eigenvalues - in the case of the third eigenvalue the difference is about 
0.01. For the higher eigenf unctions, for which kcru is sufficiently large, eg. (12. 2 p is a valid 
approximation (at kcru) without collinear resummation, we may approximate the difference, 
6uJi2, in eigenvalue between the two models in terms of the difference of the running couplings, 
oisikcriti) — Ois{kcTit2)i at the value of kcru for each model (where the frequency vanishes), 
namely 

5uJi2 ^ (asikcriti) - as{kcrit2)) — Xo(0) + ( ^'^'^"''^^^^'^ ^ 5xi(0), (4.5) 

TT \ n J 

which gives numerical results in agreement with those seen in Fig. Hlfor n > 10. 

In Fig. |5]we show a representative subset of eigenfunctions in the Standard Model and 
the SUSY model in the transverse momentum region relevant for a fit to HERA data. The 
eigenfunctions are computed at the same value of = to display SUSY effects only (in the fit 
the eigenfunctions with the same n have in general somewhat different ?7's). As expected, the 
first two eigenfunctions are identical since their values of kcru are below the SUSY threshold. 
The third and higher eigenfunctions display significant differences which affect the quality 
of the fits to data. Remarkably, these differences diminish for higher eigenfunctions and for 
n = 41 the two eigenfunctions almost overlap in this region. The reason for this can be seen 
from Fig. [2], which shows that for the relatively low transverse momenta the differences in 
the frequencies between the two models decreases with increasing eigenvalue number, so that 
if the infrared phases are equal, the functions will be almost identical in this region. 



5 Results 

One of the main results of the previous paper [T] was that we found a simple power depen- 
dence between the infrared phase f]n{ko) and the eigenfunction number n. In this paper we 
use the same functional dependence, defined for rin{ko), where ko denotes an infrared scale 
at which the phase of the leading eigenfunction vanishes. The relation between the r] phase 
at kQ and k^ is given by eq. (l3.ip . where the value of ko should be close to Aqcd- Thus we 
take the parameterization 

Vnik) = Vo , (5.1) 

V [iT-max - 1) / 

where Umax is the number of eigenfunctions we use for the fit and rjo represents the total 
range (in units of vr) of infrared phases that are used for the fit. The value of the parameters 
ko, K and rjo are determined in the fit. 
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Figure 5: Comparison of a representative subset of eigenf unctions in the Standard Model 
(blue) and the SUSY model (red) computed with the same value of r] = . The eigenvalue 
number is given in the upper right corner. 
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As explained in [T], since the eigenvalue tends to zero for large n, the form of the phase 
given by eq.f l5.ip means that as a function of eigenvalue, u, the phase has a cut singularity 
at a; = 0, i.e. 

r]{uj) = (^-) +b + cu + du'^^ (5.2) 

This allows the generalization of eq. fl5.ip by treating all the constants a, b,c,d, - ■ ■ in eq. fl5.2p 
as free parameters. We have tested these parameterizations but find no improvement in the 
quality of the fit although we introduced more parameters; we therefore used the simple 
version of the phase condition, eq.f l5.ip . in all of our fits. In ref. [1] we fixed the value of the 
parameter tjq, which represents a total range of the t] variation. In the present evaluation, 
since we have to perform many more fits, we prefer to treat it as a free parameter, to assure 
a bias free evaluation of all cases. Therefore we use in the fits the 3 parameters of eq. fl5.ip 
and the 2 parameters from the proton impact factor. For the impact factor we take the 
parameterization 

%{k) = APe-""', (5.3) 

as in ref. [I]. The fits were performed using the HERA data [2\ with x < 0.01 and > 
8 GeV^ or Q2 > 4 QeVl 

As in [Ij we find that there is an significant improvement of /Ndj in the > 8 GeV^ 
region due to various higher order effects, such as the NLO contribution to the photon impact 
factor and possibly also the proximity of the saturation region. In this region we have a total 
of 108 data points and a total of 5 parameters - so the number of degrees of freedom is 
Ndf = 103. We therefore consider the > 8 GeV^ region as our main investigation region 
and use the > 4 GeV^ as a cross check. 

We investigated the fit quality as a function of the maximal number of eigenf unctions, 
Umax- In contrast to the result of the analysis described in [1], we found that in the su- 
persymmetric analysis the best fit can be obtained with a somewhat smaller number of 
eigenfunctions; only 100 (rather than 120) eigenf unctions are required to obtain the best x^- 
For the Standard Model the best fit is obtained with rimax = 120, but only with a small 
difference in the fit quality, = 122.5(100) and x^ = 120.1(120). 

In Table [T] we show our fits for various SUSY thresholds as well as the Standard Model. 
Let us first note that the fco values obtained in the unbiased fit, ko ~ 275 MeV, are close to 
Aqcd- At the same time the value of b implies that the proton impact factor peaks around 
Aqcd, as it should be in the self consistent description. This together with the relatively 
low x^'s of all fits confirms the success of our construction of the infrared boundary. 

The quality of the fits shows a clear preference of the evaluation with SUSY effects; the 
fit for the Standard model is worse than the fits with SUSY thresholds larger than 3 TeV. 
A SUSY threshold of 3 TeV, which is close to the reach of LHC also gives a worse fit. On 
the other hand for a SUSY threshold in the region of 10 - 15 TeV, the quality of the fit is 
the best, but that for significantly larger SUSY thresholds the fit quality worsens again. 

Although the overal quality of the fit for all data with > 4 GeV^ is significantly 
worse, for reasons outlined above, the preference for N=l SUSY with the threshold region of 
10-15 TeV is also seen from a fit to the > 4 GeV^ data. In this region there are 128 
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Table 1: Fits for N=l SUSY at different scales. The bottom row corresponds to the Standard 
Model. All fits are performed with Umax = 100. 



-0.4 



-0.8 



T I r 

Standard Model 
SUSY at lOTeV 




0.25 



Figure 6: The eigenvalues and infrared phases for the Standard Model and N=l SUSY at 
10 TeV, as determined at A; = /cq- 



points and the 's of the best fits are 184.3 (3TeV), 164.5 (6TeV), 155.6 (lOTeV), 152.6 
(15TeV), 169.7 (30TeV), 164.7 (50TeV), 164.3 (90TeV). The best for Standard Model is 
169.7. The values of the fit parameters are similar to the values shown in Table [H for the 
> 8 GeV^ region. 

The consistency of the fit results and a clear preference of the SUSY fits (with a scale 
above 3 TeV) over that for the Standard Model indicates that supersymmetry improves the 
data description and suggests that some new physics similar to = 1 SUSY is present in 
the 10 - 15 TeV region. 
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6 Infrared Boundary: rj — uu relation 



The infrared boundary condition that leads to a discrete spectrum can be expressed as the 
ansatz that the phase at some infrared transverse momentum, ko, is a fixed function, rj{uj) of 
the eigenvalue u, of the form given by eq.f l5.2p . imposed by the infrared properties of QCD. 
The discrete pomeron spectrum is driven by this function. Therefore, in Fig. E]we show the 
values of the eigenvalues Un and the infrared phases ?7„ (in units of tt) both for the Standard 
Model and = 1 SUSY at a threshold of 10 TeV. The numerical values of the parameters 
on the RHS of eg. (15.21) turn out to be substantially different in the two cases. This function 
constitutes the infrared boundary conditions on the eigenvalues of the BFKL kernel. As 
explained above the eigenvalues are somewhat larger for the SUSY model, but both of these 
functions have a cut at u = 0, the order of the cut singularity being a little less in the case 
of the SUSY model (the dip is not so steep). 

As was discussed in ref. pLj the appearance of the singularity in the t] — u relation indicates 
that some important contribution to the perturbative expansion at very large transverse 
momenta is missing. Therefore it is interesting to observe that the introduction of SUSY 
softens somewhat the observed singularity (k is reduced from 0.6 to 0.56) and at the same 
time reduces the number of eigenfunctions required - it is a step towards a description of 
data using only few discrete pomerons. This could also indicate that there exist other, even 
stronger symmetries at very high energies, which are missing in the present evaluation, and 
which are responsible for the remaining singularity of the phase r]{u). 

7 Summary 

In our previous paper [1] we have shown that DPS gives a very good description of the 
low-x HERA data and it was suggested that this fit may have sensitivity to BSM physics. 
This proposed sensitivity emerged from the fact that the higher eigenfunctions have support 
over a very large range (see eq. (12.51) ). extending from the infrared region to way above the 
threshold for any new physics and that through the required phase-matching process, the 
low energy behaviour of these eigenfunctions depends on their high-energy behaviour. 

In this paper, we have shown that this is indeed the case. The introduction of N=l SUSY 
at some threshold alters the value of the /3-function and hence the rate of the running of 
the coupling. Furthermore there are contributions to the NLO characteristic function of the 
BFKL equation from the SUSY partners. Since the properties of the discrete pomeron are 
determined from a combination of the running coupling and the characteristic function, the 
eigenfunctions of the BFKL kernel are significantly affected by the introduction of SUSY. 
Notwithstanding the fact that the proposed SUSY scale is considerably above the scale 
probed at HERA, the altered high-energy behaviour of the eigenfunctions feeds into the 
low-energy, as well as generating a somewhat different spectrum of eigenvalues. 

The discrete spectrum depends on the treatment of the infrared boundary condition, 
which is now more involved due to the fact that we are using the two-loop running, 
instead of one-loop, as in [I]. Constructing this boundary we took the most conservative 
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approach of using perturbative QCD as a guideline at every step. Our previous paper 
was devoted to the task of finding the relation between the eigenfunction number and the 
phase of the eigenfunction oscillations which is essentially of the non-perturbative origin. 
Notwithstanding the substantial differences between the eigenfunctions with and without 
SUSY, we find here that the best fit is obtained using the same form of this dependence as 
was used in that paper, although other forms for this dependence were attempted without 
improving the fit quality. 

Together, the different spectrum of eigenvalues and the different shapes of the eigenfunc- 
tions in turn affect the parameters of the fit to data and also the quality of the fit. The main 
result of this paper is that if the SUSY threshold is introduced at 10 -15 TeV, the fit was 
significantly better than that of the Standard Model (x^ = HO as compared to = 122 for 
108 points and 103 degrees of freedom, i.e x^/^df = 1-06 vs 1.18). On the other hand, the 
introduction of SUSY at the threshold of 3 TeV, just within the reach of LHC, generates a 
fit which is no better than the fit obtained from the Standard Model. 

It is pertinent to emphasize the qualitative difference between the fit obtained here and 
the usual DGLAP fit [TU]. Over the low- a; region, this DGLAP fit obtains a slightly better 
value of per degree of freedom [y^ /N^f ~ 0.95). However, the DGLAP parameterization 
is designed to cover the entire range of x, whereas ours is only valid for sufficiently low-x 
where an expansion in ln(l/x) is valid. The improved quality of the fit of [lOJ is likely to 
be due to the terms with positive powers of x that are present in that fit. The important 
qualitative difference between the two fitting procedures is that the parameters obtained in 
the DGLAP approach are unaffected by any new physics at high-energy thresholds - their 
prediction for the structure functions would remain unchanged until the threshold (in Q"^) 
for new physics were reached. On the other hand, as we have emphasized in this paper, the 
values of our parameters are affected by new physics thresholds and consequently the 
evolution above the fit region will always be affected by such new physics. This considerably 
stronger predictive power of the BFKL equation is not only due to the fact that it is a 
different type of evolution equation, but also that it describes the dynamics of the gluon- 
gluon interaction which (after accounting for the infrared boundary conditions) produces the 
two-gluon quasi-bound states with a non trivial spectrum of singularities in the j-plane. 

The method described in this paper opens a new possibility to use high precision experi- 
ment to search for new physics at energy scales considerably larger than the scales at which 
the experiments are performed. It can be applied to any low-x process which was measured 
with comparable accuracy to the HERA F2 data, like the Drell-Yan or W and Z production 
at LHC. The application of this method to LHC data could lead to higher sensitivity due to 
substantially higher scales involved. The discrete pomeron solution provides a unique tool 
for such an investigation owing to the fact that the construction of the eigenfunctions is 
based on the quantum mechanical approach in which the extremely high energy (up to the 
Planck scale) and low energy behaviour of its wave functions are intimately connected. 
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